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Abstract 

We present some basic properties of two distinguished discretiza- 
tions of elliptic operators: the self-adjoint 5-point and 7-point schemes 
on a two dimensional lattice. We first show that they allow to solve 
Dirichlet boundary value problems; then we present their Darboux 
transformations. Finally we construct their Lelieuvre formulas and 
we show that, at the level of the normal vector and in full analogy 
with their continuous counterparts, the self-adjoint 5-point scheme 
characterizes a two dimensional quadrilateral lattice (a lattice whose 
elementary quadrilaterals are planar), while the self-adjoint 7-point 
scheme characterizes a generic 2D lattice. 
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1 Introduction 



In the last two decades of the XIX century and in the beginning of the XX 
century many great mathematicians (Bianchi, Darboux and others) devel- 
oped a differential geometry studying transformations of certain geometric 
structures and proved theorems of permutability of such transformations, 
obtaining in turn nonlinear superposition principles for the nonlinear differ- 
ential equations characterizing the above geometries. These results can be 
viewed as the pre-history PJ |2] of the modern theory of integrable nonlin- 
ear systems, which is based on the existence of linear differential operators 
possessing symmetry transformations of Darboux type (the Darboux Trans- 
formations (DTs)) and nonlinear isospectral symmetries (the integrable non- 
linear systems). 

After more than one century, in studying discrete integrable systems 
(which are, in a sense, richer and more fundamental than their continuous 
counterparts and, for these reasons, worth studying), one often makes use of 
the mutual interplay between geometry and the theory of integrable systems 
also in the discrete case [3]. 

The main goal of this paper is to present two examples of such an inter- 
play; we start with the self-adjoint 7-point operator 



for which the existence of Darboux transformations has been recently estab- 
lished [4j, and we obtain their geometric interpretation through the construc- 
tion of their Lelieuvre formulas. It is worth mentioning that, very often, one 
follows the opposite direction: a "geometric insight" allows one to construct 
an integrable system. 

In the above equations Tm and T„ are the translation operators with 
respect to the discrete variables (m, n) G I? 



and with the self-adjoint 5-point operator 




(2) 





and fm,n = /(^, '^) is a function of (m, n). 



2 



The classical Lelieuvre formulas [3] 

R,u = N,u X iV, R,,= Nx iV,„ . (3) 

allow one to construct a two-dimensional surface R{u, v) in from its nor- 
mal (not necessarily the unit one) image N{u,v). The co-ordinate net {u,v) 
of the surface obtained in this way is the asymptotic one and the normal field 
satisfies the Moutard equation 

N,^,= F{u,v)N. (4) 

The Moutard equation is covariant under the Moutard transformation. A 
second Darboux transformation appeared in literature (the first one was 
the transformation of Ribaucour); it gave Guichard [7j the possibility to 
describe Weingarten rectilinear congruences in a very elegant way and, in 
turn, it allowed Bianchi and other Geometers to construct many systems of 
nonlinear differential equations which now are called soliton systems. 

Another example of Lelieuvre type formula was obtained by Bianchi jH] 

R,, = N,y xN, R,y=Nx N,, . (5) 

Now the normal field satisfies the 2D Schrodinger equation 

N,.. +N,yy = FN, (6) 

which is also covariant under a Darboux transformation, and the coordinate 
net (x, y) is, in this case, the isothermally-conjugate one. 

An extension of the Lelieuvre formulas to an arbitrary co-ordinate system 
(or, better to say, to a co-ordinate free language) and to hyper-surfaces in the 
equi-affine space of arbitrary dimension has been obtained in ^ 02] • Another 
extension of the Lelieuvre formulas can be found in JT]. A discretization of 
the Lelieuvre formulas Q was first given in the paper J2] and a discretization 
of the notion of Weingarten congruences was proposed in ^2] • 

The main result of this paper consists in the construction of the Lelieuvre 
formulas for a 2D quadrilateral lattice (a lattice whose elementary quadri- 
laterals are planar) and for an essentially arbitrary 2D lattice in 
space. This result allows one to establish that the operators (j2I) and (0) 
characterize, at the level of the normal vectors, respectively, the above 2D 
lattices. 
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For our purposes it is not necessary to deal with the Euchdean space, but 
it is enough to enrich the affine space with the volume form Vol (by Vol* we 
denote the dual form of Vol); i.e., it is enough to deal with the equi-affine 
space eA'^. This enables one to construct the cross product from an ordered 
pair of linearly independent vector fields (say (a, &)); i.e., to construct the 
element N G T*eA^ such that < N\a >= =< N\b > and < N\c >= 
Vol{a] b] c} for every c G TeA^. 

The second goal of this paper consists of the illustration of some of the 
basic criteria for constructing the proper discretizations of partial differential 
operators. Again we use, as illustrative examples, the operators (H)) and ©. 

The paper is organized as follows. In section |21 we present some of the ba- 
sic criteria which we use as a guide for constructing the proper discretizations 
of partial differential operators. These criteria are systematically applied, in 
the remaining sections, to the illustrative examples given by the operators 
£5 and £7. In section 01 we show that the operators £5 and £7 preserve the 
elliptic character of their differential counterpart, being applicable to solve 
the Dirichlet problem on a 2D lattice. In section 0] we show that the opera- 
tors £5 and £7 possess, like their differential counterparts, DTs. In sections 
ini and m we derive the Lelieuvre formulas for, respectively, the continuous 
counterpart of £7, for £7, for the continuous counterpart of £5 and for £5, 
verifying that the geometric meaning of the operators £5 and £7 is the proper 
discretization of the geometric meaning of their differential counterparts. 

We conclude this introductory section with some general remarks on the 
operators £5 and £7. The operator £7 can be interpreted as the most general 
self-adjoint operator on the star of a regular triangular lattice |I6i El; it 
possesses a class of Laplace transformations [1^1 HZ] and plays a relevant role 
in a recently developed discrete complex function theory ^Hl- Its natural 
continuous limit jlj: 

Adl + + 2Sd,dy + (A,, +S,y )d, + (5„ +^,, )dy " F (7) 

is the most general, second order, linear, self-adjoint operator of the second 
order. The operator £5 is instead the most general self-adjoint operator on 
the star of a square lattice j3] and its natural continuous limit is the following 
self-adjoint elliptic (if AB > 0) operator 

Adl + A,, + Bdl + ' ^- 
It is interesting to remark that the following distinguished gauge equivalent 
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form of the operator £5 [T: 



F F F F 

P ^ m,n rp I ^ rn-l,nrp-l , m,n rp , ^ ni,n—l rp-1 /n\ 
i-SrhJnt '■— — JmH 7; J m + " "i n ~ Qm,n, ) 



r, 



m+l,n 



m,7i+l 



admits DTs and reduces, in the continuous hmit, to the celebrated Schrodinger 
operator in the plane 

dl + d'^-Q. (10) 

Therefore it can be considered as a distinguished integrable discretization of 
the Schrodinger operator 





(m,n+l) 




(m-l,n) 


(m,n) 


(m+l,n) 










(m,n-l) 





(m-l,n+l) (m,n+l) 



(m- 1 ,n) {ni,n) (m+ 1 ,n) 



{m,ii-l) (m+l,n-l 



Fig.l 



The 5 - and 7 - point schemes on the square lattice 



2 Basic criteria for discretizing partial differ- 
ential operators 

In order to construct the proper discretization of a partial differential oper- 
ator we are guided by the following criteria. 

1. It should possess a large class of (discrete, continuous, isospectral, non- 
isospectral,..) symmetries, (at least) as large as that of its differential 
counterpart. 

2. Its spectral properties should be similar to those of its differential coun- 
terpart. 

3. The discretization should preserve the hyperbolic or elliptic character 
of the partial differential operator; in particular, if the operator is ellip- 
tic, the discretization should be applicable to solve a generic Dirichlet 
boundary value problem on a 2D lattice. 
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4. If the continuous operator is geometrically significant, the discretization 
should possess a geometric meaning which generalizes naturally that of 
the continuous operator. 

In this paper we show that the difference operators (0), (j21) are discretiza- 
tions, respectively, of the partial differential operators ([7j), (jH)) which satisfy 
the properties 1., 3. and 4. (the spectral properties will be discussed else- 
where) . 

3 The Dirichlet boundary value problem 

As we pointed out in Section |2l a proper discretization of a second order ellip- 
tic operator should be applicable to solve Dirichlet boundary value problems 
on a 2D lattice. 

Consider, for the sake of concreteness, the following Dirichlet problem on 
a bounded domain of for the operator (jHJ: 



which appears very frequently in applications. 

It is easy to convince one-self that the 5-point self-adjoint scheme for the 
operator £5 (see e.q. Q): 



which, in the natural continuous limit, reduces to the above equation 
is perfectly adequate to solve a generic Dirichlet boundary value problems 
on a 2D lattice. The reasoning behind it (well-known to numerical analysts 
|19j ) is clarified by the illustrative example of Figure 2. 



{A^,, +(5^,y ),y = F^, {x, y)eVc M^ 
\E'(x,?/) given on dT>, 



(11) 




m,n—l 




o 



o 



o 



o 



o 



o 



(0,0) 



o 



o 
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Fig.2 



A simple Dirichlet problem for the 5-point scheme 



Suppose we want to solve the Dirichlet problem associated with the 5- 
point scheme (jl2|) in the subset of 7? consisting of the white and black points 
in Fig.2. If the field '^m.n is given at the boundary points (the white points), 
the unknown values of i\}ra.n at the 4 interior points (the black points) are 
uniquely constructed solving the following linear, inhomogeneous, deter- 
mined system of 4 equations for 4 unknowns: 



obtained applying 4 times the 5-point scheme ()12|) with center at the interior 
points. 

The same argument holds for more general subsets of Z^; its only possible 
failure is associated with the non generic situation in which the relevant 
matrix determinant of the system (which depends on the coefficients a, 6, /) 
is zero. 

The definitions of interior and boundary points used in the above illustra- 
tive example are intuitive: the (nearest) neighbourhood of a point (m, n) of 
the square lattice consists of the four points (m + l,?T,), (m, n + 1), (m — l,n), 
(m, n — 1). Given a subset f2 of Z^, its interior points are the points of 
for which all neighbouring points belong to its boundary points are 
instead the points of such that some of the neighbouring points do not 
belong to f2. 

We remark that the 5-point scheme (jl2j) is, among all possible difference 
equations adequate to solve Dirichlet problems on 2D lattices, the simplest 
possible scheme. 

Using similar considerations, one can show that the 7-point scheme 



is applicable to solve Dirichlet problems on a 2D lattice. Notice that, on 
a square lattice, two white points should be added to the boundary with 
respect to 5-point scheme. 



-/o.oV'o.o + «o,oV'i,o + ^0,0^0,1 = -a-i,oV'-i,o - &o -i^o -i 
ao,o^o,o - /i.oV'i.o + &i,o^i,i = -ai,o^2,o - 
^'o,oV'o,o - /coV'o,! + ao,i^i,i = -a_x,\^-\,\ - &o,iV'o,2 
^i.oV'i.o + ao,iV^o,i - /i, 1^1,1 = -ai,iV^2,i - 1^^1,2 



(13) 




(14) 
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o o o 



o ■ ■ o 

(0,1) 

o ■ ■ o 

(0,0) (1,0) 



o o o 

Fig. 3 A Dirichlet problem for the 7-point scheme on the square lattice 

4 Darboux transformations 

Non isospectral symmetries of Darboux type for linear differential operators 
play a relevant role in the theory of nonlinear integrable systems. They allow, 
for instance, to construct solutions of these nonlinear systems from simpler 
solutions, through an iterative procedure. As we mentioned in Section El a 
good discretization of a partial differential operator should preserve this type 
of symmetries. 

In this section we present the DTs for the operators £5 and Cj. These 
results are extracted from 

4.1 DTs for £5 

Consider the operator £5 together with the associated difference equation: 

where am,n, bm,n and fm,n are given functions. 

The operator £5 exhibits the following covariance property (gauge invari- 
ance) : 



•^5 ^ ^5 5'm,n-^55'n 
^m,,n ^ Oiffi fi Oim,n9m,n9m+l,ni ^m,?i ^ ^m,?i b'm,n9im,n9im,n+li (■^^) 

fm,n fm,n fm,n9m,n 

and possesses the following DTs. 

Let 6 be another solution of (fT3|) . i.e: 
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then: 



Eliminating from ()15p and p7|) we get 



{P'm—l,n'4^m,n^m—l,n (^m—l,n^m,n^m—l,n)~^ 
^n(fim,n—l^m,n^m,n—l ^m,n—l^m,ri^m,n—l) 0; 



(19) 



where 



^mfm,n ■ fm+l,n fm,n ^nfm,n ■ fm,n+l fm,n 

^—mfm,n • fm—l,n fm,n ^—nfm,n ■ fm,n—l frr 

It means that there exists a function a such that 

ntx Llm—l,n"m,n"m—l,n'-^—m a ; 



(20) 

A ^ 1, Q Q A i'm,n 

where: 

^mfm,n fm+l,n fm,n ^nfm,n fm,n+l fm,n 

^—mfm,n fm—l,n fm,n ^—nfm,n fm,n~l fm,n- 

Setting 

Tm,n n ' 

we find that tp'^ „ satisfies the following equation 



where 



Orn. ri - ' 



a. 



m—l,n 



b 9 ^'---^ = a e ^22) 

"m— l,n— 1 t'm— l,n— 1 (^m— l,n— 1 Cm— l,ri— 1 

and 

/m,n ~ ^m,,n ( Ojrn,n 7 I" ^m-l,n 7 I" ^m,?i 7 I" ^m,n-l 7\ 

\ "m+l,n "m— l,n "m,n+l "m,n—l 

(23) 

Comparing equations (fTHj) and ()23|) . we also infer that = 1/0 is a solution 

of (EID. 
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4.2 DTs for £7 

The construction of DTs presented in the previous sub-section apphes to the 
self-adjoint 7-point scheme associated with £7: 

^m,n'4^m+l,n ~l~ '^m— l,n^m— l,n ~l~ ^m,n^m,n+l ~l~ ^m,n— l^m,n— (24) 

which is a discretization of the most general second order, self-adjoint, linear, 
differential equation in two independent variables. 
Let 6m,n be another solution of equation (j21j): 

a _|_ a _ {■ a l^^J 

Eliminating from (j^^ and we get 

Ar„ ft ft (i>m_ji_ V'm-l.n N I „ ft ft ( ^rn,n~l V'm-l.n \1 i 

m[^m—l,n"m,n"m—l,n\ a a , j ~r i'm,n"m— l,n"m,n— 1 1 a , a , /J~r 

A f/l ft ft ( '4>m,n _ i>m,n-l \ I „ Q Q ( ^m-l.n _ -^m,.!-! \1 _ n. 

'-^n\'~'m,n—l"m,n"m,n—l \ a a , J ~r <'m,n"m—l,n"'m,n — l\q , a , Jl 

It means that there exists a function ip' such that 



(26) 



aft ft A 

•^m,n"m—l,n"m,n—l'-^—m a 

and function ip'^ „ satisfies the following equation 

n' 1// -I- ?// -I- ft' 1// -I- h' 7// -I- 

""771,71 ^711+1, 71 ' ""m—l,n Ym—l,n ' "m,n Ym,n+1 ' "m,n—l Ym,n—1' 

q' ih' 4- q' ?// = f 7// 

"^171+1,11 Tm+l,n—l ' "^TTtjTi+l r m— l,n+l J m,nT 171,111 

where the new fields are given by 

/ 0m,n^m-\- 1 jnC-m — 1 ,n L/ ^m,nG'm,n-\- 1 ^m,n— ■ 1 



(27) 



(28) 



/ ^m — l.Ti — — l,n"m,ri — 1 



f = ft ( n' 1 _|_ ^1 1 _|_ IJ 1 _|_ y 1 I 

i77i,7i '^"^."'^""i,ne™+i,„ ^ ^™.'^e^,n+i ^™."-ie^,n-i"^ 

■*77i+l,7ie^+l „_i ^ *m,7i+l6l„_i,„+i^ 



(29) 



10 



and where Pm,n ^m.n'^m— l,n^m,n— l ~l~ ^va—l,n^m,n(^m—l,n ~l~ "5m,n^m,n— l^m,?i— 1 • 

Again 6*^^ = is a solution of (j2Hl)- 



5 Lelieuvre formulas 

Lelieuvre idea of describing the surface parametrized with asymptotic co- 
ordinates via its CO- normal image 01201 V^^Y^ important role in the theory 
of surfaces in equi-affine spaces. Due to the generalizations of Lelieuvre 
formulas to a co-ordinate free language [3^1, one can describe the hyper- 
surface via its co-normal image in an arbitrary coordinate system. 

In this section we show that the co-normal image of a general 2D surface 
m is a vector solution of the partial differential equation 

(A*,, + {S^,y )„ +(5^,. ).y = F^ (30) 

associated with the operator ((Tj). 

We recall some basic facts. We denote by R the position vector : — >■ 
eA^ of a parametrized surface in eA? and we assume that the surface 

i) be regular, i.e.: 

N oc R,^xR,yy^O] (31) 

ii) be twice differentiable; i.e., in particular: 

Rixy Riyx ) (32) 

iii) be locally strongly convex, hence: 

Vor{N-N,^-N,y)^Q. (33) 

Then, from (|33p . we infer that the vector fields N x iV,^ and N x iV,^ are 
linearly independent and are tangent field to the surface. Therefore one can 
decompose the fields R,x and R,y as follows: 

R,y =ANx N,, +PN X N,y , 

R,x=-QNxN,,-BNxN,y, ^ ^ 

where, since R,x xR,y = {AB — PQ)V ol* {N ; N,^ ; N,y )N, due to assumption 
dSD), we have AB - PQ ^ 0. The equality < iV|^,^3^ >= gives 
(P - Q)Vor{N- N,x ■,N,y) = 0; so we have P = Q =: S and, finally, 

R,y= AN X N,x+SN X N,y, 

R,,= -SNxN,x-BNxN,y, ^ ^ 
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AB-S^^ 0. (36) 

The compatibihty condition R,xy = R,yx of equations (jH^ leads to the partial 
differential equation 

[AN,, + {SN,y +{BN,y ),y + {S N , , ),y=FN (37) 

associated with the operator (|7j), which is nothing but the most general self- 
adjoint equation of second order in two independent variables. 

Conversely, let Ni, N2 and be three linearly independent solutions of 
the self-adjoint equation which we assume to be non parabolic; i.e.: 

AB-S^^ 0. (38) 

Select any frame in eA^ and the vector field N = [Ni , N2 , N^] with respect 
to the CO- frame. Since A^i, N2 and N3 are linearly independent, we have that 
Vol*{N; N,^; N,y) 7^ 0. In addition satisfies equation (jHTj) : the vector 
multiplication of both sides of this equation for by iV it-self yields, after 
manipulation, the equation 

{AN X iV,^ +SN X N,y ),^ +{SN x N,^ +BN x N,y ),y=0, 

from which we infer that there exists a vector field R such that equations 
fl35|) hold. Interpreting R as the position vector of a surface, we infer that A^ 
is a co-normal field to this surface and that this surface is regular, since 

xR,y = {AB - S^)Vol*iN; iV„ ; N,y )N. 

6 Lelieuvre formulas associated with the 7- 
point scheme 

In the previous section we have shown that the operator (|7j) characterizes the 
co-normal image of a generic surface in eA^. According to the last criterion of 
Section 121 a good discretization of ((Tj) should possess an analogous geometric 
meaning. Indeed in this section we will show that the difference operator £7 
describes the co- normal image of a generic 2D lattice in eA^. 

Consider a lattice D U dQ — >■ eA^ and denote by fm,n the position 
vector with respect to a frame. By "lower" triangles we mean the triangles 
with vertices (rVi,„, rm+i.n, "^m.n+i) and by "upper" triangles we mean the 
triangles with vertices (r^+i,„+i, r^+i,™, r^,n+i)- 
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Fig. 4 Upper and lower triangles of the 2D lattice 

We assume that: 

• The upper and lower triangles are not degenerate; i.e., the three points 
of each triangle are not collinear. For the lower triangles of the 2D 
lattice this condition means that: 

^mXm,n ^ ^nfm,n 7^ 0. (39) 

Then we denote by ^ any co-normal non- vanishing field to the lower 
triangles: 

f)^ — A f X A r 

where is a non- vanishing scalar field. Analogously, the non degen- 
eracy condition 

Amrrn,n+1 X Anfm+l,n 7^ (40) 

for the upper triangles of the 2D lattice allows one to define any co- 
normal non- vanishing field „ to the upper triangles by: 

where is a non-vanishing scalar field. 

• The fields h^ ,^ and „ satisfy the following conditions: 

V^n ■■= yol*{n^ra,n. ^^+l,n, ^n+l) 7^ 0, (41) 

:= Vol*{hl^^, hl_,^^, <„_i) ^ 0. (42) 
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/^From assumption (jUj) it follows that the discrete vector fields x 



'^m-i,n „ X n^_^ are linearly independent and therefore they span 



the tangent space to the lower triangle (the same is true for fields n^'^m n-i 
and ^ ^m+i n-i)- So we can write 

^n'^m,n = ''^m,™ ^ i^'m-l,n'fi'rn-l,n "I" P"i~l,n^m-l,n+l ) ) (43) 

The equality < A^A„f„_„ >=< n^_„| A„A„f„^„ > is equivalent to 
{p — q)V^^^ = 0; so we have (taking into account (jH))) p = q =: s and, as a 
result, 



^m''"'m,n ^m,n ^ l^m,n— 1 ~l~ '^m,n-l'^m+l,n— l) 5 

where the coefficients a, 6, s are defined by: 

_ <"m,n+ll^"^"^ + l,"> 1 _ <"m + l.nl^"^m.n + l> 

> 

5r; 



(44) 



(45) 



^From AmA„fv« „ = A„A,mfVi,n we finally get that the lower co-normal vector 
satisfies the self-adjoint 7-point scheme 



Sm-l,nf^m-l,n+l ~^ '^ni,n-l^m+l,n-l ~ fm,n'^m,n- 

Consider now the co-vector fields: 



(46) 



Ym,n • (^m—l,nf^rn—l n ~^ '^»^— l,"-^??!— 1 n+l' 
^m,n ■ ^m,n—l^rn,n—l ~^ ^'ni,n—l'^rn+l,n—li 

then equation (j46j) can be re-written in terms of them: 
A direct calculation shows that 
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from which we infer that 



\/or(ri^,„;F^,„;Z^,„)7^0. (49) 
For the normal to the upper triangle: 



^m,n ^rn,n^n'^rn+l,n ^ ^m'^m,n+l (50) 



SO 



^ m,n ~ ■^m,n ^ m,n'^m-l,n ^ m-l,n'^m,n-l ^m,n-l ^ \"'m,m -' iri,n, ^m,n)'^ 
iflm,n^m—l,nbrn,rJ^rn,n—l ~l~ ^m,n—l^m—l,n'^m,rJ^rn~l,n ~l~ ^m,n~lbm—l,n'^m,n^m,n—l)' 

(52) 

Therefore 

0'm,n^m—l,nbm,n^m,n—l ~l~ ^rn,n—l^m—l,n^rn,n^m—l,n ~l~ ^m,n~l^m~l,n^m,n^m,n—l 7^ 0- 

(53) 

Summarizing, we have the following theorem. 

Theorem 1. Consider a two-dimensional lattice Z 9 f2 — >■ eA^ such that its 
position vector frn,n and its lower fi^ n ^'^^ upper co-normals obey the 
conditions / T^) . ^Jl^ and J^^ . T/ien there exist functions am,n, bm,n 

and Sm,n obeying conditions / f37]) and / f33|) swc/i t/iat i/ie following Lelieuvre 
type relations hold 



^n'^m,n ^m,n ^ ('^m-l,"^m-l,n "I" '^m--l,n.^m-l,n+l J ' (54) 



and swc/i t/ifl^ i/ie lower co-normal field satisfies the 7-point self-adjoint scheme 

(^m,nflrn+l,n ~^ ^'ni-l,n'flm-l,n ~^ ^m.,n'^?Tt,n+l ~^ ^rn,n-lft'ra,n-l'^ (55) 
Sm-l,nf^rn-l,n+l ~^ '^ni,"-l'^m+l,n-l ~ fm^rfi'm^n- 

Conversely, consider the field ri^„ satisfying: i) equation \5^) with the 
coefficients obeying conditions \51]) and \5'J^) : ii) the conditions and 
\49^ . Then the Lelieuvre type formulas \54^ define the position vector rm,n 
of a 2D lattice in cA^, having as a lower co-normal. The position 

vector, the lower co-normal and the upper co-normal given by ■ — 
A^,,Amfm,„+i X A„fm+i,„ Satisfy the conditions (E^j, and j^^. 



m,n 
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7 Lelieuvre formulas associated with the self- 
adjoint 5-point scheme 

In the previous two sections we have shown that, on the level of the Lelieu- 
vre type description, the operator ((Tj) and its discretization £7 characterize 
respectively a generic 2D co-ordinate net and a generic 2D lattice in eA^. In 
this section we introduce distinguished reductions on the above generic nets 
(lattices), showing that i) the reduction from a generic net to a conjugate 
net (a surface parametrized by conjugate co-ordinates) is characterized, on 
the level of the Lelieuvre type description, by the reduction from the general 
self-adjoint partial differential operator ((Tj) to the self-adjoint operator (jH)); 
ii) the reduction from a generic 2D lattice to a 2D quadrilateral lattice (a 
lattice whose elementary quadrilaterals are planar) is characterized, on the 
level of the Lelieuvre type description, by the reduction from the 7-point 
scheme (I14j) to the 5-point scheme (jl2j) . 

Let R{x, y) be the position vector of a conjugate net and let iV(x, y) be 
any co-normal vector field; then: 

< N\R^ >= =< N\Ry >, definition of N, 

< N\R^y >= 0, conjugacy. 

Therefore, assuming that the surface be locally strongly convex (condition 
(j33|) ). one infers that the following Lelieuvre type formulas hold: 

R^ = BNy X N, Ry = AN X N^, (56) 

where 

B(.,y)= . A(.,y) = ^mL^. (57) 

Vol'(N,N,,N,) Vol-(N,N,,N,) 

The integrability condition R^y = Ry^ implies equation 

{AN,, + A,N, + BNyy + ByNy) X iV = , (58) 

which is equivalent to the wanted self-adjoint equation 

{AN,), + {BNy)y = FN, F = F{x, y). (59) 
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Conversely, it is straightforward to prove that, if N satisfies equation (|59p for 
some coefficients A, B, F, then the vector R, defined by (j^ . is the position 
vector of a conjugate net having N as normal vector. 

For the discrete case we follow the same reasoning. Any co-normal vector 
field of a quadrilateral lattice satisfies equations: 

< Nm,n\^mrm,n >= =< iV,„,„ | A„fm,„ >, definition of N, 

< Nm,n\^m^nrm,n >= 0, quadrilaterality. 
Therefore, assuming that 

Vo/*(iV^,„, iV„+i,,, iV^,„+i) ^ 0, (60) 

the following Lelieuvre type relations exist between the tangent vectors and 
the co-normal to the lattice 

^m^m,n ^m,n— l-^m,n ^ -^m,n—ly ^nTm,n (^m—l,n-^m,n ^ -^m—l,ni 

(61) 

where the scalar fields am,n and bm,n are defined by: 

> 



Vol''(Nm,n,Nm + l,n,N^^ri + l)'' ('gO') 
l,n l^m'^mjri + l 



'm,n Vol*{Nm,n,Nm+l,„,Nm,n + l)' 

The integrability condition for equations (jF)T|l implies the equation 

('^m,nA^m+l,n + '^m-l,nA'm-l,n + ^m,nA^m,n+l + &m,ri-l-^m,n— l) X iVm,n = 0, (63) 

which is equivalent to equation C^N = 0. 

Conversely, it is also easy to show that, if Nm,n satisfies equation (jHSl), 
then the Lelieuvre formulas (pT|) define a proper embedding of a 2D quadri- 
lateral lattice having n as co-normal. 
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Fig. 5 The 5-pomt scheme for the normal vector 

We remark that this result could have been deduced in a faster way, f.i., 
from ^^Et ). observing that the reduction from £7 to £5, expressed by the 
equation Sm,n = 0, is equivalent to the condition that the tangent vectors 
of the upper triangles of the 2D lattice are perpendicular to i.e., it is 

equivalent to the condition that the 2D lattice is quadrilateral. 

We conclude this section remarking that, due to the gauge covariance 
properties of the operators ^ and £5 (see (IT^ and its continuous limit), the 
normal vectors appearing in the characterizing equations and C^N = 
have an arbitrary normalization. The characterizations (see, e.g., pT] ) 

N^y + + ^Ny = 0, 

of respectively a 2D conjugate net and of a 2D quadrilateral lattice in terms 
of their co-normal vectors are instead gauge dependent. 
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